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Abstract. The Cauchy problem for harmonic maps from Minkowski space with its stan- 
dard flat metric to a certain non-constant curvature Lorentzian 2-metric is studied. The 
target manifold is distinguished by the fact that the Euler-Lagrange equation for the ener- 
gy functional is Darboux integrable. The time evolution of the Cauchy data is reduced to 
an ordinary differential equation of Lie type associated to SL(2) acting on a manifold of 
dimension 4. This is further reduced to the simplest Lie system: the Riccati equation. Lie 
reduction permits explicit representation formulas for various initial value problems. Addi- 
tionally, a concise (hyperbolic) Weierstrass-type representation formula is derived. Finally, 
a number of open problems are framed. 
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1 Introduction 

Let (M, g) and (N, h) be Riemannian or pseudo-Riemannian manifolds and (p : M — > N a smooth 
map. The energy of ip over a compact domain PCMis 



where A is the Laplacian on M and K 1 ^ the Christoffel symbols on N. A map <p is said to be 
harmonic if it is critical for e((p). Harmonic maps generalise harmonic functions and geodesies 
and have been under intense study since the pioneering work of Eells and Sampson [8]; see [3] 
and [12] and references therein for comprehensive introductions to the field. If domain(<^) = 
M = M then harmonic maps are geodesic flows. If codomain(<^) = iV = M. then harmonic maps 
are harmonic functions. If (M, g) is pseudo-Riemannian then harmonic maps are called wave 
maps. In this paper our focus is on wave maps, specifically, the case (M,g) = (R 1 ' , dxdy). Any 
further reference to wave maps in this paper means the domain space M is Minkowski space Mr' 1 



This paper is dedicated to Peter Olver on the occasion of his 
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especially his influential, diverse applications of Lie theory. 




The critical points of e((p) satisfy the partial differential equation (PDE) 
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with its standard flat metric. To provide slightly more insight, a "physical" illustration of wave 
maps in this class can be given: the motion of a frictionless elastic string constrained to vibrate on 
Riemannian surface (N, h), such as a sphere, is exactly described by a wave map into TV; see [17]. 

There is a well-known geometric literature on wave maps that has developed over the last 
two decades, especially their existence as solutions of completely integrable systems; see [11] for 
a textbook account with many references. There is a closely related physics literature where 
the relevant systems are known as nonlinear sigma models; see [22]. The first person to treat 
the Cauchy problem for wave maps into Riemannian targets was Chao-Hao Gu [10]. He estab- 
lished the fundamental result that for smooth initial data, wave maps into complete Riemannian 
metrics have long-time existence. Gu's work initiated many further investigations where regu- 
larity constraints on the initial data have been significantly relaxed. Furthermore some higher 
dimensional problems have been treated; see [16]. 

In this paper we initiate the study of the Cauchy problem for wave maps in the special case 
where the systems they satisfy are Darboux integrable. Our first main result proves that the 
solution of the Cauchy problem for such a Darboux integrable nonlinear sigma model can be 
quite explicitly expressed as the flow of a special vector field £&i,fc 2 (see Theorem 1) which itself 
is a curve in a certain Lie algebra of vector fields canonically and intrinsically associated to 
the Darboux integrable nonlinear sigma model, namely, its Vessiot algebra. In consequence of 
this, standard constructions which facilitate the resolution of systems of Lie type such as Lie 
reduction become available to the solution of the Cauchy problem for such wave map systems. 
For this reason we have included an appendix to this paper which gives a brief summary of the 
main results on systems of Lie type adapted to the applications we envisage. In this paper we 
have decided to focus on just one interesting nonlinear sigma model in order to discuss the rela- 
tionship between Darboux integrable hyperbolic systems on the one hand and the resolution of 
the corresponding Cauchy problem via differential systems of Lie type and to do so as explicitly 
as possible. However, it will be seen that the proof of Theorem 1 is easy to generalise to other 
Darboux integrable systems. Indeed the very recent work [1] outlines a general, intrinsic proof 
of the close relationship between systems of Lie type and the Cauchy problem for a wide class of 
Darboux integrable exterior differential systems. The second main result of the paper marshalls 
the general theory of Darboux integrable exterior differential systems [2], and generalised Gour- 
sat normal form [18, 19] to derive a hyperbolic Weierstrass-type representation (Theorem 3) for 
wave maps into the non-constant curvature metric (1). 

As the name implies the notion of Darboux integrability originated in the 19th century and 
was most significantly developed by Goursat [9]. Classically, it was a method for constructing 
the "general solution" of second order PDE in one dependent and two independent variables 

F{x, y, U, U x , Uy, U X xi U X yi Hyy) — 

that generalised the so called "method of Monge" . It relies on the notion of characteristics and 
their first integrals. We refer the reader to [9, 13, 20, 21] for further information on classical 
Darboux integrability. There are also extensive studies of Darboux integrable systems relevant 
to the equation class under study in the works [14] and [23]. 

In this paper we use a new geometric formulation of Darboux integrable exterior differential 
systems [2]. At the heart of this theory is the fundamental notion of a Vessiot group which, 
together with systems of Lie type are our main tools for the study of the Cauchy problem for 
wave maps. 

The PDE that govern wave maps have the semilinear form 

u xy = / (x, y, u, u x ,u y ), u,f£ R n . 

Each solution possesses a double foliation of curves called characteristics. Such PDE often 
model wave-like phenomena and projection of these curves into the independent variable space 
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describes the space-time history of the wave propagation. The characteristics of u xy — f are 
the integral curves of a pair of rank n + 1 distributions 

H X = {D X + D y f ■ d Uyy , d Uxx } , H 2 = {D y + D x f ■ d Uxx , d Uyy } , 

where 

D x = d x + u x d u + u xx d Ux + fd Uy , D y = d y + u y d u + fd Ux + u yy d Uy , 

are the total differential operators along solutions of the PDE. Note that if is the standard 
Cartan codistribution for u xy = f then H\ © H 2 = ann#. Distributions Hi are well-defined 
with canonical structure. 

Definition 1. If A is a distribution on manifold M then a function / : M — > IR is said to be 
a first integral of A if Xf = for all X E A. 

There is a geometric definition of Darboux integrable exterior differential system [2]. For 
wave map equations it reduces to 

Definition 2. A semilinear system u xy = f with u, f € M. n is Darboux integrable at a given 
order if each of its characteristic systems Hi has at least n + 1 independent first integrals at 
that order. 



2 The Cauchy problem 

We consider wave maps 



u 



(Ri-^,,/ w )^(R» (1) 



R. Ream [15] studied the PDE for wave maps into nonzero curvature surface metrics that are 
Darboux integrable on the 2- and 3-jets and proved a theorem that any such metric is (real) 
equivalent to one or other of the metrics 

du\ + du\ 
1 ± e ui 

Here we consider a semi-Riemannian version of a Ream metric and study the corresponding 
Cauchy problem. We show how the solution of the Cauchy problem for wave maps (1) can be 
expressed as an ordinary differential equation of Lie type. Indeed we prove that the solution of 
the Cauchy problem for wave maps is naturally equivalent to an initial value problem of a Lie 
system for a local action of SL(2) on a manifold that is locally diffeomorphic to M 4 . This is further 
reduced to an initial value problem for a single Riccati equation together with a quadrature. 

The target metric in (1) does not have constant curvature nevertheless is globally defined and 
positively curved everywhere; in fact K = 2 _1 (1 + e Ul )~ l . However, the metric is nonetheless 
very special because the wave map system turns out to be Darboux integrable, as demonstrated 
in [15]. 

The Lagrangian density for this metric is 

_ Ul x Ul y ~ U 2x U 2 y 

~ 1 + e-"i 
whose Euler-Lagrange equation is 



u\ x ui y + u 2x u 2y ui x u 2y + u 2x u ly 

2(l + e«i) " ' U2xy 2(l + e"i) 
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The change of variables {u\,U2) i-> ((ui + it2)/2, (ui — U2)/2) = (u,v) transforms this to 

'Ui'x, '^y ^y / \ 

u *v + 2 (x + eU /2+«/2) = °' ^ + 2 (1 + e «/2+,/2) = °- ( 2 ) 

We now prove 
Theorem 1. Consider the initial value problem 

u *v + 2 (l + e «/2+«/2) = °' ^ + 2 (1 + e «/2+*/2) = °' W 

where j is a curve with tangents nowhere parallel to the x- or y-axes, n is a unit normal vector 
field along 7 and fa, ipi are smooth functions along 7. 

1. Problem (3) has a unique smooth local solution. Moreover, the unique local solution is 
expressible as the solution of an ordinary differential equation of Lie type associated an 
action o/SL(2) on M 4 . 

2. Given the unique local solution (u, v) from part 1, the Cauchy problem for harmonic maps 

is given by 

U\ = u + v, U2 = u — v. 
where u\, U2 satisfy initial conditions 

du du 
ui\ y = fa + 02, u 2 \ y = fa - fa, -q^- =ipi+ip2, ~Qn 1 = ^ 1- ^ 2 - 
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Proof. By hyperbolicity, the problem is locally well-posed. System (2) has four first integrals 
on each characteristic system, H\, H2. Let us label the first integrals 

y, (3i, fa, /S3 for Hi, and x, ct\, a 2 , 03 for H 2 . 

For this system a%, f5\ are first order differential functions while 012, 03, /?2, are of second 
order. Finally, while the 8 first integrals are functionally independent, we have 

da\ d(3i 
—— = a 2 and —— = p 2 - 
ax dy 

Let k\(y), fc 2 (y) be arbitrary functions and consider the overdetermined PDE system defined 
by (2) together with the additional PDE 

Pi = k!{y), fo = k{y), fo = Wv), (5) 

where the dot denotes ^-differentiation. It can be shown that this overdetermined system £' 
is involutive and moreover admits a 1-dimensional Cauchy distribution, as we will see. Now 
suppose we fix a smooth curve 7 embedded in a portion of the xy-plane, M and suppose Cauchy 
data is prescribed along 7 as in the Theorem statement. Then by an argument similar to [21, 
Proposition 3.3], 7 can be lifted to a unique curve 7 : Kl-> J 1 (7V r ,IR 2 ) which agrees with 
all the Cauchy data. Let l : %kiM ~ > J 2 (J\f,M 2 ) denote the submanifold in J 2 (J\f, M 2 ) defined 
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by PDE system £' and the contact system on J 2 . Let 0} tl k 2 = t*S be the Pfaffian system 
whose integral submanifolds are the solutions of £' . Our aim is to choose the functions k±, k2, 
if possible, in order that we can extend 7 to a 1-dimensional integral submanifold 7 of #fci,fc 2 - 
Now 

0fci,fc 2 = W 1 = du — u x dx — Uydy, 0J 2 = dv — v x dx — v y dy, 00 3 = du x — u xx dx — ^dy, 

oo 4 = du y — ^dx — Uyydy, cj 5 = dv x — v xx dx — f 2 dy, a; 6 = dv y — f 2 dx — v yy dy}. 

Pulling back by 7 we observe that a; , co 2 pullback to zero by construction. Forms ui and w 6 
define the functions u yy and v yy along 7 while to 3 , w 5 define functions u xx and v xx along 7. 
All these functions are expressed in terms of the Cauchy data, (pi, tpi. Substituting these back 
into (5) determines the functions £7, ^2 in terms of (pi, For completeness we give the first 
integrals of H2 

_ u x v x _ dai 

011 ~ 1 + exp ( - S±*) ' ° 2 ~ dx ' 

(2f x u xa; - 2u x v xx - u x v\ + v x ul + 2(f x u xa; - it x v xx - u x v 2 x + ti 2 ^) exp (^)) 
013 ~ u x v x (l + exp(^)) ' 

Those of Hi are similar but with y replacing x. 

We will now use these first integrals to demonstrate that Oki,k 2 has a one-dimensional Cauchy 
distribution and that in particular the Cauchy vector can be chosen to be a curve in a certain 
Lie algebra - the Vessiot algebra [2] of system (2). It will be seen that the Cauchy vector is 
generically transverse to the Cauchy data and extends the one-dimensional integral submanifold 
of Ok ± .k 2 t° the solution of the Cauchy problem. Because the Cauchy vector is a curve in a Lie 
algebra, this extension from a one-dimensional to a two-dimensional integral of ^fc li fc 2 is an 
ordinary differential equation £ of Lie type. Its coefficients and initial conditions are fixed by 
all the data present in the problem, including the Cauchy data. Any solution of £ (independently 
of its initial conditions) permits a Lie reduction of £ and will permit us to solve the IVP for £. 

Indeed, setting z\ = e"/ 2 , Z2 = e"/ 2 , Z3 = u y , Z4 = v x , a,i = c^, bi = fy, i = 1, 2, 3 we calculate 
that the contact system on J 2 (M 2 ,M 2 ) pulled back to PDE (2) is ^ = {k 1 , . . . , k 6 }, where 

1 a a 2 . . 3 a 01(1 + 31*2), ^iz 3 

k = da\ — a<idx, k = db\ — 02dy, k = dz\ dx dy, 

ZZ2Z4 2 

4 , Z 2Z4, b 1 (l + Z 1 Z 2 ) 5 dlZ-i 1 2 > 

k =dz 2 —dx dy, k = dz 3 + dx + -(z 3 - b 3 z 3 - b\)dy, 

2 2z\z 3 ZZ1Z2Z4 2 

61 1/2 \ r b\Z/± 

k = dzi H — (za — (I3Z4 — ai)dx H dy. 

2 2z 1 z 2 z 3 

Pulling back to submanifold (5) yields a Pfaffian system with 1-dimensional Cauchy distri- 
bution spanned by 

=d y - k ^{R 1 + AR A ) + ^1r 2 + ^3, 

where the Ri form a basis for the Vessiot algebra 1 

Ri = z\d Zx - z 2 d Z2 - 2z 3 d Z3 , R 2 = 1 + ZlZ2 d z + d Z3 Z4 —d Z4 , 

ZlZ 3 z 1 z 2 z 3 

R 3 = z\z 3 d zi - z 2 d Z3 , R A = -- (ztd Zl - z 2 d Z2 ) 



1 A brief geometric construction and interpretation of the Vessiot algebra is given in Section 4 and Appendix A. 
See [2] for a complete exposition. However, detailed knowledge of Vessiot algebras is not a prerequisite for this 
paper. 
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with nonzero structure 

[ill, R2] = 2i?2j [Rl, R3] = — 2i?3) [R2,R&\=Rl- 

Since Ckik2 * s a curve in the Vessiot algebra it determines an ODE of Lie type. Furthermore 
£ki,ka is generically transverse to the Cauchy data. 

Note that p\ = R\ + 4i?4, p 2 = R2, P3 = R3 generates a local action of SL(2) on M 4 : 

[pl,p 2 ] = 2p 2 , \pl,P3] = -2p3, [P2,P3] = Pi 

and the Cauchy vector is 

ZkiM = d y~ k 2{y)pi + h(y)p2 + 2^3- 

In summary, vector field £,ki,k 2 flows the 1-dimensional initial data solution curve 7 of 0^ fc 2 to 
a 2-dimensional solution. This completes the proof of Theorem 1. ■ 

Example 1. As an illustrative example we consider the initial value problem 



UxUy V x Vy 

2(1 + e «/2+«/2) = U ' v *y + 2(l + e«/2+«/2) 
<9it <9t> 



mi = ui = 0, — — = -— = y/2, 

where 7 = Since x and y are light-cone coordinates x = (£ + r)/2, y = ((, — t)/2, the 

curve 7 corresponds to time r = 0. Thus we have constant initial values at time r = 0. We wish 
to determine the system of Lie type whose solutions corresponds to the solution of this Cauchy 
problem. We have n = 2~ 1 / 2 (d x — d y ) and we get 

U X \^ = V x{j = 1, U y^=Vy^=-l. 

So our initial curve in J 1 is j(x) = (z, x, 0, 0, 1, — 1, 1, — 1) = (x,y,u,v,u x ,u y ,v x ,v y ). This 
translates to an initial curve in the adapted coordinates 

(x, y, zi, z 2 , z 3 , z A ) = (x, x, 1, 1, -1, 1). 

We extend this to a unique 1-dimensional integral of \& and get 

7(2;) = (x,y,z 1 ,z 2 ,z 3 ,z 4 ,ai,a 2 ,a 3 ,b 1 ,b 2 ,b 3 ) = ^x, x, 1, 1, -1, 1, i, 0, 0, ^, 0, 0^ . 

Thus we get 

fc 1 (y) = l/2, fca(y) = 0. 
The Cauchy vector is therefore 

^,0 = ^ + ^2 + ^3. 
We flow this vector field obtaining solutions Zi(x, y) subject to the intial conditions 
zi(x,x) = 1, z 2 (x,x) = 1, z 3 (x,x) = —1, Z4(x,x) = 1. 
The ODE are 

dzi 3^2 1 + ziz 2 dz 3 1 Z4 

ZlZ 2 , 



dy ' <9y 4ziz 3 4 2 ' <9y 4ziz 2 ^3 
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to be solved for Zi(x,y) subject to the given initial conditions. In fact, we need not solve the 
whole system but only the equation for z 3 and then substitute this into the equation for z\ . This 
gives the function u up to a quadrature after which the function v can be obtained algebraically 
from the PDE itself 

e v/2 = _ e -u/2 I 1 



2u X y 



We get the unique solution of the Cauchy problem to be 

( y/2 V2 V 

u(x, y) = v(x, y) = In ( cosh ~^( x — y) + V% sinh -^-(x — y) j . 

Example 2. A slightly more interesting example is obtained from the initial conditions 

, ^ 1 1 dv 

n, 7 =21nA, v u =2ln-, = = 1, 

for any A > 0. The Cauchy vector is £1 and its flow subject to the initial conditions along 
y = x being 

Zl = \ z 2 = \, * = -js> ** = ^ 
gives rise to the unique solution 

+ y\ exp (-|) (3 + 2V2) - exp (-§ ) \ 



x 

u = 2 In I A exp ' 



2(1 + v^) 



n1 (l (x + y\ exp(-y) (2^2 + 3) -exp (-a) (2^-3) -2exp(-^) 
v = 2 In — exp 1 1 



2A *A 4 J exp(-|) {1 + V2) +exp(-f) (l - V§) J' 

Thus, even constant initial data has the potential of producing some interesting explicit solutions. 
Indeed, one can ask if this solution is global in time. 

The fact that we only had to solve for z 3 in Example 1 (and Example 2) holds not only for 
these illustrative examples since the system of Lie type in the general case is 

dzi 1 dz 2 . , , (1 + z x zz) 

- -Z1Z3, -t— = ki(y)- 



dy 2 dy 2z x z 3 

■q- = o( fc i (y) + k 2 {y )z 3 -z 3 ), — = -- . 

dy 2 V o; dy 2z 1 z 2 z 3 

This proves 

Theorem 2. In the Cauchy problem for wave maps 



(R^dxdy) - (M,^=^t) 



0/ Theorem 1, let functions k\(y) and k 2 {y) be defined as described above. Denote by £ i/ie initial 
value of Z3 = u y along 7. Let T be the unique function satisfying Riccati initial value problem 

r)T 1 

- = -(k 1 (y) + k 2 (y)T-T*), T {y = (. (6) 
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Uxy + — 7~uT7T~77rr = 



Define function u to be the unique solution of 

du , 
^ = r ' ^ = ^ 
and let v be defined algebraically from the partial differential equation 

UxUy 

2 (1 + e «/ 2 +«/2) 

upon substituting solution u and solving for v. The functions u, v constitute the unique solution 
of the Cauchy problem. 

This implies that the solution of any given Cauchy problem for wave maps into the metric 

<iu? — du% 

h = — 

1 + e" - "! 

relies on the solution of a Riccati initial value problem together with one quadrature. The 
interesting point here is that the Riccati equation is the simplest non-trivial equation of Lie 
type. It is a Lie system for the local SL(2)-action on the real line that globalises on MP 1 . 
In general, solutions of Riccati equations develop singularities in finite time, even those with 
constant coefficients. However, the theorem above provides a correspondence between Cauchy 
data for the wave map and the Riccati initial value problem (6). An interesting problem is to 
study this correspondence and link the nature of the Cauchy data with the properties of the 
solution of (6) and in turn, link this correspondence with the geometry of the target metric. 

For the standard initial value problem where Cauchy data is posed along = 2r = x — y, the 
relationship between the Cauchy data and the coefficients of the Riccati equation is complicated. 
However, due to its significance and for latter use, we give it explicitly: 



1 fa (y/2fa x ~ 2fa) exp ( f J ! (V2fa x + 2fa){y[2fa x + fa) exp ( f 

h = --r 7 — t~t\ — > °i 



1 exp + exp (f ) ' 4 exp + exp (- ^ 



2 



k 2 = (T 1 ^2 [fa* + 2ip 2 2 - 2V2il> 1( f> lx - 4feiJ (^ 2 + spl& 

+ 4^01 (4> lx - fa y/tj exp (-^ - ^) + 8 (~fa x ^2 + (fa + V2fa x 

where 

5 = A(V2fa x + fa) (fa x -faV2) . 

Recall that functions fa, fa are respectively the values of u and v along 7, the fa are the values 
of the normal derivatives along 7 as stated in the theorem. For instance for arbitrary constant 
initial conditions along the time axis 

fa(x) = u(x, x) = 2 In Ai, ^2(2^) = v(x,x) = 2 In A2, 

fa = aV2, fa = pV2; X u A 2 > 0, 

we have 

a/3\i\ 2 , , {a-(3)\ 1 \ 2 . , 

kl{v) = T+A^' HV) = ~ 1 + A X A 2 • (7) 

The solution of the Cauchy problem in this case is easily calculated as we did in Examples 1 
and 2 because it amounts to solving Riccati equation (6) with constant coefficients ki(y), k 2 (y) 
given by (7); indeed Example 1 is the choice Ai = \ 2 = 1, a = (3 = \[2. However the general 
formula is complicated and of itself not very informative so we refrain from recording it here. 
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3 Cauchy problem for wave maps and Lie reduction 

We have shown that to solve the Cauchy problem for the wave map equation it is enough to 
"evolve" the initial data curve 7 by solving a Riccati equation which is notably the simplest 
nontrivial system of Lie type. One significant feature of Lie systems is that they admit "reduc- 
tion by particular solutions", otherwise known as Lie reduction. There are only a few sources 
scattered in the literature on this topic, among them [4, 5, 7]. In this section we give an illus- 
tration of how Lie reduction may be useful in resolving instances of the Cauchy problem for our 
wave map system. Appendix A. 2 of this paper summarises the known results on systems of Lie 
type, oriented toward the applications at hand. In this and subsequent sections we will adopt 
the notation and theory set out in Appendix A. 2, to which we refer the reader. 
Consider the Cauchy problem for wave maps with (non-constant) Cauchy data 

01 = 02 = 0, Vi = -2V2, V2 = 2V2x 



along the curve y = x. We get 



ki{y) = -ty, k 2 {y) = -y-l, 
so that the corresponding Riccati initial value problem of Theorem 2 is 
BV 1 

— = -(-2y-(y + l)T-T 2 ), T(x,x) = 2. (8) 

We observe that r = 1 — y is a solution vanishing at y = 1. Implementing the procedure 
described in Appendix A. 2 obtains one factor in the fundamental solution 



9o(y) 



1 1-y 
1 



A curve in the isotropy group of is denoted by H and has the form 

The curve of Lie algebra elements associated to the Riccati equation 

dz 

— = ct (y) + 2ai(y)z - a 2 (y)z 2 
ay 



is 



For the Riccati initial value problem (8) we have 
A{y)=(~T y+l) 

V 2 4 / 

and the reduced fundamental equation is 



dy 
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where (see Appendix A. 2, especially Theorem 6 &; Appendix A. 2.1) 
B{y)=go{y)- 1 A{y)g {y)-g,{y)- ld9 ^=(^ y rl U 



dy v \ -b + fy' 

valued in the isotropy subalgebra at 0, as expected. The ODE initial value problem for the 
fundamental solution factor gi(y) is 

^-/3i7i = 0, ^ -^ 7l + /3 l72 = , 71(1) = 1, 72 (1) = 0, 

where Pi = {y — 3) /4. This problem can be explicitly solved in terms of elementary functions 
giving 



7i (y) = exp 



(y-l)(y-5) 



72(2/) = v^exp (^y - y 2 - y^j ferfi(l) + erfi Q (y - 3)^ , 

where i = y— 1 and erfi denotes a concomitant of the error function: 
2 /"^ 

erfi(x) = — ^= / exp (i 2 ) di. 

This data enables us to construct the fundamental solution g(y) = go(y)gi(y) for (8) and leads 
to the solution of the ODE in (8) 

r(y;g) = \( y )(l), 

where Ah denotes the linear fractional transformation (12) by h £ SL(2). 

It is now a simple matter to determine the value of q that satisfies the initial condition 
T(x,x) = 2 and giving the unique solution T(x, y) of the Riccati initial value problem (8). We 
find 

p, x 1 2exp(|y-ly 2 ) 
T(x,y) = l-y 



A(x)-v^exp(-|)erfi(l(y-3))' 



where 



, 2exp(-|a; + lx 2 ) ^ ( 9\ p ,(\, 
A(x) = ^ x) 4 } + v^exp (--j erfi (-(x - 3) 

Finally, we obtain an integral representation of the solution u satisfying 
du 



dy T(x,y), u(x,x) = 0, 



namely 



ry 

u(x, y) = I T(x, s) ds 



or in terms of spacetime coordinates (£, r) 

■i«-r)_ fi 



u(£,t) = / r -(f + T),s (fa. 
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The significance of Lie reduction in our ability to solve the Riccati equation should here be em- 
phasised. Without this, it would have been impossible to construct the fundamental solution and 
there would be no hope of constructing function T and constructing the integral representation 
of the solution of a Cauchy problem with non-constant initial data. 

In the example above we relied on knowledge of a simple solution, namely T(x, y) = 1 — y 
to perform the reduction. But even with polynomial or rational coefficients a Riccati equation 
will not generally have any rational solutions. In this case however we can appeal to the well 
known fact that every Riccati equation can be linearised. 

Lemma 1. The general Riccati equation 
dz 

— = ao(t) + 2ai (i)z(t) - a 2 (t)z(tf (9) 
can be transformed to the form 

±y{r) = P{r) + y{r)\ 
where z = p(t)y(t) 

Inp = J 2a>i(s)ds, r = J at2(s)p(s)ds. 

Provided these quadratures can be carried out then the explicit solvability of (9) depends on the 
properties of its linearisation 

d 2 ip 1 d 

^2 + P(r)il> = 0, where y(r) = ^V'(t)- 

Any solution of the 2nd order linear ODE can be used in the Lie reduction of the Riccati equation. 

As a consequence of Lemma 1 and Theorem 7 of Appendix A. 2, we have 
Corollary 1. The solution of the Cauchy problem for wave maps of Theorem 1 
.1,1 v ( du l-du 2 2 



(M 1 ' 1 , dxdy) -> \ M, 



1 + e~ Ul 



is reducible to quadrature provided a particular solution of the Riccati equation (6) is known. 
A particular solution of (6) can be constructed by quadrature and the solution of a linear second 
order ODE. 

Remark 1. As a consequence of Lemma 1 and Theorem 7 the differential equations solver in 
MAPLE - dsolve is very often able to construct an explicit representation in terms of known 
special functions to a Riccati initial value problem when the coefficients are polynomial functions 
of the independent variable. 



4 Hyperbolic Weierstrass representation 

In this section we use the Darboux integrability of the wave map equation (2) to compute its 
general solution and hence construct a hyperbolic Weierstrass-type representation for wave maps 
into the corresponding metric. According to [2], we pull back f to a suitable integral manifold 
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Mi, M2 of h[°°^ and i^ 00 -* respectively. It is convenient to define Mi by y = b\ = b 2 = ^3 = 
and M2 by x = a\ = a 2 = a 3 = 0. This gives Pfaffian systems 

\I>1 = I dzi — - (a 3 z 4 + a\ — zT) dx, dz 3 + - — a±z 3 dx, dz 2 - \z 2 z 4 dx, 

[2 ' 2z\z 2 z 4 2 



dz\ ai(l + z\Z2\dx. da\ — a 2 dx . 

2z 2 Z 4 ' J 



1 



and 



^2 = { dzi - \z 3 z\dy, dz 2 - 61 (1 + ziz 2 )dy, dz A + - — hz 4 dy, 

2 2ziz 3 2z\z 2 z 3 



dz3 - ^(hz3 - zl + b\)dy, dh\ - b 2 dy\, 



each of rank 5 on 8-manifolds: (Mj, Locally Mi has coordinates x, z\, Z\, z 2 , z 3 , 24, a\, a 2 , 
a 3 while M2 has local coordinates y, z\, z\, z 2 , z 3 , 24, b\, b 2 , 63. Using these local formulas, we 
define a local product structure 

(m x M 2 ,§i 0$ 2 ) , 

where 

$1 = { dq A - - (a 3 q 4 + a\ - qj) dx, dq 3 + — - a±q 3 dx, dq 2 - -q 2 q 4 dx, 

I 2 2q\q 2 q 4 2 

dqi ai(l + q\q 2 )dx, da\ — a 2 dx 

2q 2 q 4 



and 



^2 = {dpi - \pip 3 dy, dp 2 - —^—hil+pip^dy, dp 4 + — L hp 4 dy, 

2 2pip 3 2pip 2 p 3 



dp 3 



^{b 3 p 3 ~P 3 + h)dy, dh - b 2 dy\. 



As described in [2], the relationship between Vf, *$>i and ^ 2 is that every integral manifold 
of ^ is a superposition of an integral manifold of ^1 and *$> 2 . The superposition formula is the 
map 2 

(T : Mi x M2 -> M 
defined by 

1 u\ I 1 ~~ PW2P4 + ^p\P2PAqz - p\P2qz - 93 + piP2)qi 

cr(x, p, a; y, q, b) = 



PiP2 

P2(2p 4 g3g2gl ~ 919293 ~ P49291 + 9l92 ~ Pi + 1) 

9391 

(2y>493 - 93 - Pa + l)p 3 p\p 2 

1 - pip 2 p 4 + 2p 1 p 2 p 4 q 3 - p\p 2 q 3 - 93 + Pip 2 ' 



2 See Appendix A for further details on the superposition formula; in particular, how it is defined and con- 
structed. 
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(2P4<73 - <73 - Pi + l)<7l<?2<74 u u u 

, x, y, ai, a 2 , a 3 , &i, 2j »3 



2P4<?3<72<71 - 919293 - ^4^291 + 9192 ~ £>4 + 1 

(zi, z 2 , z 3 , z 4 , x, y, ai, a 2 , a 3 , 61, 6 2 , 63) ■ 



The usefulness of this factorisation of the integration problem for is not only that the inte- 
gration of \Pj relies on ODE while that of relies on PDE but that the are locally equivalent 
to prolongations of the contact system on J 1 (M, M 2 ). To see this we turn to the characterisation 
of partial prolongations of such contact systems provided by [18, 19] which also provide simple 
procedures for finding the equivalence. To implement this we compute the annihilators 

ann := H± 

la , a , a l0- + I1Q2) a . <?2<?4~ «ig3 o . !/ , 2\ o a q 1 

= < o x +a 2 d ai ^ d qi + ——d q2 - d q:j + -(a 3 g 4 + a x - q 4 )d q4 , d a2 , d a3 }, 

[ 2(/2<?4 2 lq\q 2 q 4 2 j 

ann ^2 := -£^2 

, , bl(l+Piy2)a 1/, 2 , i\a &1 ^ 4 » a a 1 

= j e y + b 2 d bl + — a pi + v 2yip3 9 P2 + 2 (&3P3 - p 3 + - ^wps ^ fe2 ' 63 J ' 

We show that each of the -ffj is locally equivalent to the partial prolongation C(0, 1, 1) of the 
contact distribution on J 1 (R,M 2 ). That is, the contact distribution on J 1 (M,IR 2 ) partially 
prolonged so that one dependent variable has order 2 and the other order 3 with canonical local 
normal form 

C(0, 1, 1) = {d t + z\d z i + z\d z x + z\d z 2 + z\d z -i + z\d z ^ d z i, dq}. 

Let T> be a smooth distribution on manifold M and assume that T> is totally regular in the sense 
that T>, all its derived bundles and all their corresponding Cauchy bundles have constant rank 
on M. Denote mi = dimPW, y* = dim Char and x\~i = dim Char T>^_ 1: where 

Char V[% = V^ l) n Char V (i) . 

Below, k denotes the derived length of D. 

According to Theorem 4.1 of [18], a totally regular distribution D on smooth manifold M is 
locally equivalent to a partial prolongation of the contact distribution on J 1 (M,IR' ? ) for some q 
if and only if 

1. The integers mj, x 3 , Xj-i satisfy the numerical constraints 

X 3 = 2m,j — rrij + i — 1, < j < k — 1, Xi-i = m» — 1, 1 < i < k — 1. 

2. If mfc — mk-i > 1 then a certain canonically associated bundle called the resolvent is 
integrable 3 . 

A pair (M, I?) that satisfies these conditions is said to be a Goursat manifold or Goursat 
bundle. Moreover, if "D is a Goursat bundle on M then it is locally equivalent to a partial 
prolongation with x 3 ~ Xj-i dependent variables at order j < k and nik — mk-i dependent 
variables at highest order, k. This uniquely identifies the partial prolongation associated to 
a given Goursat manifold. Before discussing the question of constructing equivalences, let us 



3 In the original formulation of Theorem 4.1 in [18], the integrability of Char 23j_j is an additional hypothesis 
to be checked. This is a simple task but unnecessary since it is easy to see that this bundle is always integrable 
and the hypothesis can be omitted. 



14 



P.J. Vassiliou 



solve the recognition problem for the distributions Hi. We take H\ which we denote temporarily 
by K. We find 

Char£ = {0}, Char = {d a2 , d a3 }, Char k[ 2) = {8 a2 , 8 a3 , d ai , d q3 }, 
Char = {d ai , d a2 , d a3 , d qs , d x }. 

Calculation shows that the dimensions of the derived bundles are 

dimK = 3, dim^ (1) =5, dim^ (2) = 7, dimK (s) =8, s > 3. 
Hence the derived length is k = 3. Below we check the first condition of a Goursat bundle. 
Table. Checking the numerical constraints satisfied by (M,T>). 



j 


rrij 


rrij-i — 1 


2rrij — rrij+i — 1 




X J 





3 




6-5-1=0 







1 


5 


2 


10- 7 - 1 = 2 


2 


2 


2 


7 


4 


14 - 8 - 1 = 5 


4 


5 



Hence, K is a Goursat bundle with k = 3, m/- — mk-i = 1 and the only nonzero difference 



X'' 



Xj-i being at order j = 2: x 



xi 



5 — 4 = 1. Hence there is one variable of order 2 and 



one variable of order 3. This solves the recognition problem and we can assert that K := H\ is 
locally equivalent to C(0, 1, 1). Next, we show how to construct an equivalence. Given a Goursat 
bundle, an efficient method for constructing an equivalence map was worked out in [19] and relies 
on the filtration induced on the cotangent bundle. Denote by H^) and Sjli the annihilators 
of Char K^) and Char K^\, respectively. Then we obtain the filtration 

S (2) c s (2) c 
spanned as 

{dqi, dq 2 , dq 3 } C {dqi,dq2,dq 3 ,dx} C {^1,^2,^3,^,^4,^1}- 

The construction proceeds by building appropriate differential operators and functions. Because 
m k — m k-i = 1, condition 2 of the definition of Goursat manifold is vacuous. Instead we fix 
any first integral of Char T>^ k ~ l \ denoted, t and any section Z of T> such that Zt = 1. Then, 
define a distribution H k inductively as follows: 



[Z, U e ], n 1 = Char V^, l<£<k-l. 



There is a function ip k which is a first integral of Tl k such that dip k A dt 7^ 0. The function ip k 
is said to be a fundamental function of order k. The space of fundamental functions of lower 
order can be constructed from the filtration above by taking quotients. Specifically, as noted 
above, in this case the only fundamental functions of less than maximal order 3 are of order 2. 
They are described by the quotient bundle 



:f) / s (2) 



{dx} mod dq\ , dq 2 , dq% . 



Without loss of generality we can take </? 2 = x to be a fundamental function of order 2. The 
construction of the equivalence map is now as follows. Function t is the "independent variable" 
and successive differentiation gives the higher order variables 



V, 



Z <p , 
Z (f , 



4 



v3,„k 
Z (f . 
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We now implement this. The first integrals of Char are spanned by qi, q2, 93 and any 

function of these can be chosen to be t. If we choose (say) t = q%, then for Z we choose 



9294 



where X is the first vector field in the basis for K above; for then Zt = 1, as required. We then 
construct the integrable distribution II 3 as described above and discover that its first integrals 
are spanned by 

A "?l93 

q2 and 



1 + 9192 

Hence the fundamental function of highest order is 

2 k 9193 

z = cp — 



1 + 9192 

The data 

, 1 2 9193 
t = q 2 , z = x, z = — 

1 + 9192 

and differentiation by Z now constructs the local equivalence ip identifying K = H\ and 
C(0, 1, 1). The local inverse ip^ 1 : K — > Mi determines the integral submanifolds of Hi. 

An exactly analogous calculation holds for H2 and one arrives thereby at an explicit map 
if>2 '■ K — > M2 representing the integral manifolds of H2- The explicit integral manifolds of 
are a superposition of those of Hi and H2 : 

RxI-> <r(V'r 1 (M),V'2" 1 ( K ))- 

In this way we obtain remarkably compact representations for wave maps into this metric: 

Theorem 3 (hyperbolic Weierstrass representation). For each collection of twice continuously 
differ entiable real valued functions fi(s), f2{s), gi(t), g2(t) of parameters s, t, the functions 

x = fi(s), y = gi(t), 

e u/2 = {tg2{t)-l)f2{s) + {2tg2(t)-l)f2{s) 2 

92(t)(f2(s) + Sf 2 (s)) 

v /2 = (gMg) - 1M) ± - ^{tf 

f2{s){92{t) + tg 2 {t)) 

define harmonic maps 

(R^^xdy) -> (n, d ^ + ~J^ by Ul =u + v, u 2 = u-v. 



5 Concluding remarks 

We've seen that for a large family of non-constant initial data, it is possible to construct explicit 
integral representations for solutions of the Cauchy problem for wave maps into a certain non- 
constant curvature metric due to the fact that the corresponding Euler-Lagrange equation is 
Darboux integrable and because the Cauchy data can be extended as a flow by a system of 
Lie type. We have also constructed a hyperbolic Weierstrass representation for such wave maps 
making use of the general theory in [2, 18, 19]. The fundamental ingredients throughout include 



16 



P.J. Vassiliou 



the theory of systems of Lie type and the notion of a Vessiot group associated to any Darboux 
integrable exterior differential system [2]. We expressed the evolution of the Cauchy data as 
a system of Lie type for the action of a subgroup of the Vessiot group. 

We may perhaps regret the occurence of integral representations in our solution of the Cauchy 
problem preferring the elimination of all quadrature. Alas, this is surely a forlorn hope in 
a difficult nonlinear problem, especially when it is recalled that even in the general Cauchy 
problem for the (l+l)-linear wave equation quadrature cannot be eliminated, according to the 
d'Alembert formula. However, one can ask if there are Darboux integrable nonlinear sigma 
models with solvable Vessiot groups for harmonic maps into nonzero curvature metrics. This 
would make the application of the theory of systems of Lie type very useful indeed. In fact there 
is at least one such sigma model [6]. 

Interesting open problems include: what intrinsic properties of a metric render the corre- 
sponding wave map system Darboux integrable? Moreover, what do we learn about the geometry 
and topology of target manifolds from so vast a reduction in the Cauchy problem? This is espe- 
cially intriging when it is recalled that the solution to the Cauchy problem has been expressed 
as a curve or "evolution" in a finite Lie group, arising by Lie's theory from a corresponding 
curve in its Lie algebra. 



A Appendix 

A.l The superposition formula 

For completeness, in this appendix we make a remark on the construction of the superposi- 
tion formula <x : M\ x M 2 — > M. A general construction valid for any decomposable exterior 
differential system was worked out in [2]. Below we present results for the wave map system 
studied in this paper. The hyperbolic structure of the wave map system in adapted coordinates 
is given by 

H = H 1 (BH 2 , 

where 

u \a , a , a i( 1 + Ziz 2 ) a , Z2Z4 aiz 3 a ,1, . 2 \a 
Hi = <d x + a 2 d ai H d Zl + -—-d Z2 - d Z3 + -(a 3 z 4 + a x - z%)d Z4 , 



9(12 1 ( • 



H 2 = ^dy + b 2 d bl + — d Zl + 2ziZ3 d Z2 + -(6 3 ^3 " Z, + b t )d Z3 ~ ^^ 4 , 

db 2 , db 3 

Calculation shows that the infinitesimal symmetries of Hi which are tangent to the level sets of 
all the first integrals x, y, a, b, the tangential characteristic symmetries of Hi, are spanned by 

£ i = \ 0(21*30*1 - z%d Z3 ), +ZlZ2 d Z2 + d Z3 ^— d Z4 , -\{zid Zl + z 2 d Z2 ) + z 3 d Z3 , 

Z ZiZj, ZiZ 2 Z^ i 

~{zid Zl + z 2 d Z2 ] 
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Similarly, the tangential characteristic symmetries of H 2 are spanned by 

(l + z x z 2 z 3 1 

£2 = < -z o Zl d Z3 + -d Z4 , 

[ 2z 2 z 4 Z1Z2Z4 2 

z 2 z A d Z2 - zjd Z4 , ^(zid zi - z 2 d Z2 ) + z 4 d Z4 , --{zid zi - z 2 d Z2 ) 

The structure of these Lie algebras are in reciprocal relation, namely, for E\ the nonzero Lie 
brackets are 

[ei,e 2 ] = e 3 , [ei,e 3 ]=ei, [e 2 ,e 3 ] = -e 2 , 

while for E 2 we have 

[/l,/ 2 ] = -/ 3j [/i,/s] = -/i, \h,h\ = h and [e i ,f j ] = ) Vi,j. 

Hence, £±, £ 2 define the isomorphism class of the Vessiot algebra for this Darboux integrable 
system, in the terminology of [2]. Each £j is isomorphic to sl(2)©IR and each frames a neighbour- 
hood of a point of M 4 . By the converse of Lie's second fundamental theorem there is a local Lie 
group (G,m), where m : G x G —> G denotes group composition such that £{ coincide with the 
infinitesimal left and right translations on G. From the expressions for either E\ or £ 2 , by com- 
puting flows or otherwise, we can compute function m. The components of m coincide precisely 
with the first 4 components of the superposition map cr, where we interpret ({pi}f = i, {qi}i = i) 
as local coordinates around (e, e) G G x G; e being the identity in G. 

Remark 2. The tangential characteristic symmetries [20, 21] described above do not in general 
lead directly to the Vessiot algebra, which is needed for the construction of the superposition 
formula. In general a sequence of coframe adaptations described in [2] is required. In terms of 
these the distribution Hi ® H 2 above can be constructed from components of the fifth adapted 
coframe of a so-called Darboux pair. However, it turns out that for relatively low dimensional 
examples like the one in this paper, the full machinery of [2] can sometimes be avoided and 
instead it is convenient to carry out a direct computation of the tangential characteristic sym- 
metries of H\ and H 2 as we did above. A simple example where this procedure is already not 
sufficient is provided by a Fermi-Pasta-Ulam equation (see [20]), uu = u~^u x 



'• X X ■ 



A. 2 Systems of Lie type 

Let (a:GxM-^Mbea left-action of a Lie group G on manifold M and Q C X(M) its Lie 
algebra of infinitesimal generators; X(M) is the Lie algebra of all smooth vector fields on M. 
Let Xi be a basis for Q. Then, a vector field of Lie type or Lie field is a curve in Q 

X = Y J oa{t)X i , 

i 

where are smooth functions of parameter t. An ODE of Lie type is the differential equation 
determined by a Lie field 

dt lx{t) 

We pause to recall the map ju : g — > X(M), defined by 
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where (f e (u) is the flow (in G) of u G and g is the matrix Lie algebra associated to G, viewed as 
a matrix group. The map p, is the standard anti-homomorphism, induced by the left-action /z, 
from the matrix Lie algebra q to the Lie algebra of smooth vector fields X(M) on M. Remark 
that for each x G M, ker£ti is equal to the subalgebra X x C Q of infinitesimal generators which 
vanish at x; that is, ker/x^ is the isotropy subalgebra at cc. An action /i is effective on M if 
the global isotopy group I(M) is trivial. Recall that X(M) is the set of all g G G such that 
fj,(g, x) = x for all x G M. It follows that if the action is effective then ker /x is trivial and 
the map fi, is injective. Henceforth we assume that the G-action is effective or at least almost 
effective; an action is almost effective, if Z(M) is discrete. In this case, fl is injective for those 
elements of g whose flows are close to the the identity element in G. 

We will also assume that the action is transitive so that for each x, y G M the corresponding 
isotropy subalgebras are isomorphic, I x — T y . We therefore speak of the isotropy subalgebra or 
subgroup. We now recall some useful notation from [7]. 

Let / : M — > G be a smooth map and let u = 5(f) = df • f . We know that a; is a right- 
invariant Maurer-Cartan form and that it is valued in the Lie algebra q of G. 

Definition 3. Let M be a smooth manifold and G a Lie group. 

1. A differential form oj G £l(M, g) is said to be a Maurer-Cartan form if it satisfies 

(iw+wAw = 0. 

2. A map / : M — > G is said to be an integral of a Maurer-Cartan form oj if and only if 
5(f) =u. 

Theorem 4. On a simply connected manifold every Maurer-Cartan form has an integral. 

We denote the set of all q- valued Maurer-Cartan forms on M by Cl(M, q). Specialising to 
the case M = Mwe have the following fundamental result. 

Theorem 5. Let fj, : G x M — >■ M be an effective and transitive left- action of Lie group G on 
manifold M. Suppose uj = A(t) dt G f2(K,g) is smooth at t = 0. The unique solution of the 
initial value problem of Lie type 

dx 

^=£W)W x(0)=qeM 
is the function x q : R — > M defined by 
x q(t) = Kd(t),q), 

where g : M — >■ G, is the fundamental solution. That is, 5(g) = uj and g(0) = e. 

Remark 3. This theorem reduces the construction of the solution of a system of Lie type to 
constructing the fundamental solution g : E — > G satisfying 

f t =A(t)g(t), g(0) = e. 

This problem can still be very challenging. Lie's approach is to simplify the problem by making 
use of any known solutions (Lie reduction). See [7] for a proof of Theorem 5. 

To explain this, continuing with transitive and effective left-action fi : G x M —> M, suppose 
that xi(t), . . . , Xk(t) are k particular solutions of the Lie equation (10), satisfying 

Xj(t ) = qj, l<j<k, 
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for some collection of points qj £ M. Denote by G q the isotropy subgroup of q G M, 



Gn 



{g e G\g ■ q = q}. 



Let 



H — G qi n G q2 n • • • n G qk , 



and let f) be the Lie algebra of H. Furthermore, let a curve go : I C 
Kgo(t),qj) ■= go(t) ■ qj = Xj{t). 



G be defined by 



(11) 



Such a curve of group elements is defined up to a multiplication on the right by a curve in the 
the joint isotropy subgroup of the initial conditions. Constructing go{t) involves the solution of 
algebraic equations with a potentially large solution space with no canonical choice of solution. 
In general one aims to find a solution go(t) which passes through the identity at parameter 
value to where initial conditions are to be posed. 

Let C°°(M, G) denote the set of smooth maps from M to G. As in [7] define gauge transfor- 
mations p(h) : tt(M,g) -)• £l(M,g) by 

P {h)u) = Ad(h)u> + 8{h), Wh e C°°(M, G). 

Theorem 6 (Lie reduction; see [4, 7]). For any oj = A(t) dt 6 f2(R, g), 

1. wi = p (goit)- 1 ) lj G 0(R,f)). 

2. // gi : / CK — )■ H is an integral of uj\ then go(t)gx(t) is an integral ofu. 
A. 2.1 Application of Lie reduction to the Riccati equation 

The Cauchy problem for the integrable wave map system has been reduced to an initial value 
problem for Riccati equation (6) together with one quadrature. In this subsection we briefly 
illustrate the use of Lie reduction in the resolution of the general Riccati initial value problem. 
It is used in Section 3 of the paper. The standard action of the special linear group G = SL(2, R) 
on the real projective line induces a left-action on the real line 



a b 
c d 



(12) 



"0 


f 




"1 " 




"0 


0" 





0_ 




-1 




1 






the linear fractional or Mobius transformations. 
With basis 

it_i 

for s[(2,M) it is easy to see that the anti-homomorphism A : sl(2,] 

\((XqU-i + a\UQ + oi2Ui) = (ao + 2ct\x — a2X 2 )d x . 
Hence 



X(M) can be expressed 



A(t) = a (t)u-i + a\(t)uQ + a 2 (t)ui 
satisfies 



ai(t) a (t) 
a 2 {t) —ai(t) 



X (A(t)) = (ao(t) + 2ai(t)x - a 2 (t)x 2 )a a . 
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That is, A(t) C sl(2,M) is associated with the Lie equation 
— = ao(t) + 2a\(t)x — a2(t)x 2 . 

Now suppose that a solution xo(t) of equation (13) is known and suppose that xo(0) = 0. Solving 
equation (11) we obtain 



(13) 



90® 



1 x {t) 
1 



The isotropy subgroup of G E is 
a 



H = G 



c a 



from which the initial value problem for the fundamental solution gi (f ) of the reduced equation 
of Lie type can be deduced to be 

dgi 
dt 

where 



B(t) 9l (t), gi(0) 



(14) 



«1 — Xo«2 
— «2 







-(«i - x a 2 ) 



dt. 



B(t)dt = 5( gi (t)) = p go(t) -i ( w ) 
Since 51(f) is a curve in H, we have 
91 (*) 

With /3(f) = ai — xoa2 5 (14) has solution, expressed in terms of quadrature 



71 (i) 
72(f) 7i(i) _1 



7i(i) = exp/ P(t)cLt, 72(f) = — 7i(*) 1 / 02(^)71 (r) z 



dr. 



Finally, the unique solution of the initial value problem x(0) = g £ 1 for the Riccati equation 
can now be expressed in terms of quadrature: 



= 9o(t)gi(t) ■ q 
That is, 



7i(*) + »o(t)72(t) xo(*)/7l(t) 
72(0 71(f)- 1 



x (t) = - 7l(t)g + x „(t). 

w 72W9 + 71(f)- 1 

Theorem 7. iyef a?o(f) &e a known solution of a Riccati equation (13). Then by Lie reduction 
the construction of the fundamental solution is reducible to quadrature. 
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